Anisotropic magnetism, superconductivity, and the phase diagram of Rbi_2;Fe2-j;Se2 
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We report the crystal growth and structural, magnetic, conductivity, and specific heat inves- 
tigations of Rbi_3;Fe2-j;Se2 single crystals with varying stoichionietry prepared by self-flux and 
Bridgman methods. The system exhibits a strongly anisotropic antiferromagnetic behavior below 
400 K. Bulk superconductivity is found in samples with Fe concentrations 1.53 < 2 — y < 1.6, 
whereas for 2 — y < 1.5 and 2 — y > 1.6 insulating and semiconducting behavior is observed, re- 
spectively. Within the measured range of variation of the Rb concentration (0.6-0.8) no correlation 
between the Rb content and the lattice parameters of the samples was found. The superconducting 
samples show the smallest value of the lattice parameter c compared to the non-superconducting 
samples. The sharpest transition to the superconducting state, the highest transition temperature 
Tc of 32.4 K, and the highest diamagnetic response corresponding to a critical current density jc 
of 1.6 X 10* A/cm^ (at 2 K) is found for compositions close to Rb2Fe4Se5. Upper critical fields 
Hc.2 of ~ 250 kOe for the in-plane and 630 kOe for the inter-plane configurations are estimated 
from resistivity studies in magnetic fields. In the non-superconducting samples with the Fe con- 
centration below 1.45 both specific heat and susceptibility revealed an anomaly at 220 K which 
is not related to antiferromagnetic or structural transformations. Comparison with the magnetic 
behavior of non- superconducting samples provides evidence for the coexistence of superconductivity 
and static antiferromagnetic order. 

PACS numbers: 74.70.Xa 74.62.Bf 74.25.Ha 74.25. Bt 



I. INTRODUCTION 

Iron-based superconductorsi"— are currently among 
the most intensively studied materials. Among different 
groups of these superconductors, the iron chalcogenides 
have recently attracted particular attention. The un- 
doped iron selenide, FeSe, exhibts a relatively low critical 
temperature Tc ~ 8 K at ambient pressure,"* but it raises 
to 37 K under external pressure<^i^ Earlier attempts to 
increase T^. of FeSe using chemical doping resulted in a Tc 
~ 14 K by substitution of Se with TeJ^ Recent reports 
of Tc ~ 30 K in potassium intercalated FeSe^ further 
pushed the research activity in the iron chalcogenide fam- 
ily. Consequently, successful intercalation of other alkali 
metals (Rb and Cs) in FeSe was realized and supercon- 
ducting (SC) samples with Tc between 27 and 33 K were 
preparedji2ri£ Further studies of SC chalcogenides with 
hypothetical stoichiometry Ao.8Fe2Se2 (A = K, Rb, Cs, 
Tl) revealed significant differences in their SC proper- 
ties compared to the related SC pnictides with a similar 
structural arrangement. 

Intriguing coexistence of superconductivity and 
static antiferromagnetic order and a proximity to 
an insulating state were suggestedJ^"— In addition, 
angle-resolved photoemission spectroscopy (ARPES) 
studies22ri2i showed a different topology of the Fermi sur- 
face in Ao.8Fe2Se2 compared to other Fe-based super- 
conductors suggesting a pairing mechanism distinct from 
s± symmetry.^'5 Despite numerous reports, the data on 
Ao.8Fe2Se2 materials with regard to their intrinsic prop- 
erties is far from being conclusive. For example, the con- 
ducting and magnetic properties of the Rb-based samples 



reported by different authors vary significantly indicating 
a strong dependence on preparation conditions and impu- 
rity content j -'^^d'^d^d^d^i^Q Moreover, the correlation be- 
tween the properties and stoichiometry of Ao.8Fe2Se2 has 
been not fully established yet. Here we present the results 
of the structural, magnetic, conductivity, and thermody- 
namic characterization of the Rb-Fe-Se system performed 
on single crystals grown by two different methods: self- 
flux and Bridgman techniques. The selection of the Rb- 
based system is motivated by an ease to grow large single 
crystals with a higher volume fraction of the SC phase 
compared to K- and Cs-based Fe chalcogenides. The vari- 
ations of the conducting and magnetic parameters with 
the stoichiometry determined by wave-length dispersive 
x-ray electron-probe microanalysis (WDS EPMA) are 
summarized in a phase diagram of the Rb-Fe-Se system. 



II. RESULTS AND DISCUSSION 

A. Preparation and composition analysis 

Polycrystalline FeSe synthesized from the high-purity 
elements (99.98 % Fe and 99.999 % Se) and 99.75 % Rb 
were used as starting material for the growth of single 
crystals. Handling of the reaction mixtures was done in 
an argon box with residual oxygen and water content less 
than 1 ppm. The preparation conditions are given in Ta- 
ble m In the growth runs with starting composition cor- 
responding to nominal stoichiometry Rbo.8Fe2Se2 the so- 
lidified ingots showed significant inhomogeneity. Samples 
from the top of the ingots revealed superconducting prop- 



erties whereas the bottom of the ingots contained poorly 
crystahized material with ferromagnetic behavior. A sim- 
ilar separation was observed for the growth runs with the 
starting stoichiometry Rbo.9Fe2Se2, although the sam- 
ples from the top of the ingot were not superconduct- 
ing. In the growth runs with the starting stoichiometry 
Rbo.8Fe1.gSe2, such a strong phase separation was not ob- 
served, but the samples were also non-superconducting. 



The deviation from the nominal stoichiometry is a 
general problem for the intercalated iron chalcogenides. 
In most earlier reports, the composition of samples was 
determined by semi-quantitative EDX analysis or more 
precise ICP analysis. Both methods of compositional 
analysis revealed a different stoichiometry compared to 
the nominal Rbo.8Fe2Se2f^>^'^ Various ranges of iron 
stoichiometry for the superconducting samples were re- 
ported, but the accuracy of these methods appears not 
to be sufficient to detect subtle changes in the composi- 
tion. Therefore, the composition of the grown samples 
was determined by WDS EPMA that allows to reach an 
accuracy of 0.5 % for Fe and 1 % for Se. The concen- 
tration of the elements in the samples was measured us- 
ing a Cameca SX50 analyzer on freshly cleaved samples 
minimally exposed to air ('^ 1 min) to prevent oxida- 
tion. The EPMA data are also given in Table HI They 
represent the values averaged over multiple (at least 10) 
measured spots with an area of 80 x 60 /im^ . The concen- 
trations of Rb and Fe elements are normalized assuming 
a Se concentration of 2 per formula unit. The EPMA 
analysis revealed notable deviations in the composition 
of the samples from the starting stoichiometry for all pre- 
pared batches. The concentration of Fe in the supercon- 
ducting samples from different batches was found to vary 
between 1.53 and 1.60, while the variations of the Rb 
concentration were between 0.65 and 0.8. The Rb de- 
ficiency is probably caused by a strong interaction with 
quartz glass at high temperatures. The stoichiometry of 
the samples exhibiting the highest SC parameters (from 
batches BR 16 and BR 26 described below) was close 
to Rbo.8Fe1.6Se2 (or Rb2Fe4Se5 in other notation). The 
concentration of Fe in the samples prepared from exactly 
the same nominal mixture Rb2Fe4Se5 was, however, less 
than ^ 1.5 and the samples were not superconducting. 
Another group of non-superconducting samples was also 
found, with the Fe concentration being higher than ~ 1.6. 
The variations of the Fe and Se concentrations on the in- 
dividual samples were less than 1 and 2 %, respectively, 
indicating their high homogeneity. The maximal varia- 
tions of the Rb concentrations were of about 5 % which 
may indicate an inhomogeneous distribution of Rb ions. 
However, we do not exclude that this can be also an arti- 
fact due to observed sensitivity of the counting rate to a 
slight misalignment of the sample due to remaining frag- 
ments of the removed layers on the surface after cleavage 
and strong sensitivity of Rb to oxidation. 
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FIG. 1. (color online) X-ray diffraction patterns for different 
superconducting Rbi_2,Fe2-i/Se2 samples. Inset: image of 
the single crystals BR 26 subtracted from the ingot. Sample 
labels BR and F stand for those grown by Bridgman and flux 
methods, respectively. 
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FIG. 2. (color online) X-ray diffraction patterns for different 
non-superconducting Rbi_a;Fe2-ySe2 samples. 



B. X-ray diffraction and structure 

Figures 1 and 2 demonstrate the room temperature 
x-ray powder diffraction patterns obtained on crushed 
single crystals using a STOE Stadi P diffractometer with 
Cu Kq radiation, A = 1.54056 A for superconducting 
and non-superconducting samples, respectively. The lat- 
tice parameters of the samples were obtained from the 
Rietveld refinement using the FULLPROF SUITE2£ as- 
suming tetragonal symmetry I4/mmm. They are pre- 
sented in Table HI The refinement was done using the 
site occupations for constituent elements as determined 
by EPMA. 

It is important to note that non-superconducting sam- 
ples from the batch BR 19 with Fe concentration ~ 1.5 
show the narrowest diffraction lines. X-ray spectra of the 
non-superconducting samples from other batches show 
additional features which are probably related to the dif- 



TABLE I. Preparation conditions, EPMA, and x-ray data for Rbi- 
deviations.) 



rFe2-j/Se2 samples. (Numbers in brackets indicate standard 



Sample 


Preparation conditions 


Concentration 


Starting mixture 


Lattice constant | 


Method 


Soaking 
temperature 


Soaking 
time 


Cooling rate 


R,b 

(l-x) 


Fe 

(2-y) 


Se 


a, b 

[A] 


[A] 


F 266 


self-flux 


1030° C 


3 h 


6°C/h 


0.653(30) 


1.532(8) 


2.000(35) 


Rb-|-2.5FeSc 


3.9285(3) 


14.6003(18) 


F 274 


self-flux 


1030° C 


3 h 


6°C/h 


0.672(15) 


1.430(8) 


2.000(17) 


Rb+2FcSe-|-0.5Se 


3.9270(9) 


14.6155(49) 


F 286 


self-flux 


1030° C 


3 h 


6°C/h 


0.609(22) 


1.608(9) 


2.000(33) 


Rb+2.5FcSe4-0.5So 


3.9180(11) 


14.6273(61) 


F 295 


self-flux 


1030° C 


3 h 


6°C/h 


0.724(18) 


1.560(10) 


2.000(19) 


Rb-f2.5FeSc 


3.9310(7) 


14.5914(36) 


BR, 16-1 


Bridgman 


1070° C 


3 h 


3 mm/h 


0.796(28) 


1.596(8) 


2.000(37) 


Rb+2.5FeSc 


3.9228(7) 


14.5909(38) 


BR 16-8 










0.801(28) 


1.602(9) 


2.000(37) 




3.9304(6) 


14.6050(28) 


BR 16-10 










0.759(42) 


1.601(9) 


2.000(47) 








BR 17-1 


Bridginan 


1070 °C 


3 h 


3 mm/h 


0.737(26) 


1.439(8) 


2.000(38) 


0.9Rb+2FeSo 


3.9117(5) 


14.6516(36) 


BR 17-4 










0.733(24) 


1.439(9) 


2.000(38) 








BR 17-5 










0.733(24) 


1.432(15) 


2.000(23) 








BR 17-10 










0.786(17) 


1.459(10) 


2.000(15) 








BR 18 


Bridgman 


1070 °C 


5 h 


3 mm/h 


0.659(32) 


1.537(9) 


2.000(36) 


Rb-|-2.5FeSc 


3.9282(4) 


14.5899(23) 


BR 19 


Bridgman 


1070° C 


5 h 


3 mm/h 


0.698(17) 


1.497(11) 


2.000(25) 


Rb+2FeSe-f0.5Se 


3.9379(1) 


14.5983(7) 


BR 22 


Bridgman 


1070° C 


5 h 


3 mm/h 


0.644(31) 


1.620(12) 


2.000(35) 


Rb+2.5FcSe+0.5So 


3.9021(7) 


14.6484(58) 


BR 26 


Bridgman 


1070° C 


5 h 


3 mm/h 


0.740(36) 


1.600(6) 


2.000(25) 


Rb-|-2.5FeSc 


3.9190(9) 


14.5515(57) 


BR 28 


Bridgman 


1070° C 


5 h 


3 mm/h 


0.685(21) 


1.588(9) 


2.000(40) 


Rb-|-2.5FeSc 


3.9251(6) 


14.5932(36) 



ference in their vacancy structure. For example, sam- 
ples with the Fe concentration below 1.45 exhibit much 
broader peaks indicating inhomogeneity and stress. Sam- 
ples with the Fe concentration above 1.61 exhibit split- 
ting of peaks which can be attributed to orthorhom- 
bic distortions. We also noticed that superconducting 
samples exhibit much broader x-ray diffraction peaks 
and additional reflexes beside the main peaks. In the 
superconducting samples from batches F 295, BR 16 
and BR 26, the additional reflexes are of a consider- 
able magnitude. Additional peaks, although of smaller 
magnitudes, were also observed in the single crystal x- 
ray diffraction spectra in a number of SC Ao.8Fei_a;Se2 
compounds but at lower angles compared to the main 
peaks F^iii^ii^i^ These additional peaks were attributed 
either to inhomogeneous distribution of the alkali ions 
or to a modulated structure resulting from ordering of 
iron vacancies. The absence of such peaks in the non- 
superconducting samples probably exclude any signifi- 
cant inhomogeneous distribution of the Rb ions in our 
samples. Recent single crystal x-ray and neutron diffrac- 
tion studies of a number of Ao.8Fe2Se2 compounds iden- 
tified their crystal structure within a tetragonal I4/m 
space group symmetry] ^^ i^^" — The x-ray patterns of our 
samples were also refined within the I4/m symmetry. 
For example, for the SC sample BR 16 the obtained 
value of the lattice constant a=8.7977 A agrees well 
with v5a of the supercell. We must note that x-ray 
pattern only of the non-superconducting sample BR 19 
with the sharpest diffraction peaks is fully compatible 
with the I4/m space group symmetry. However, within 
the I4/m symmetry it was not possible to refine ad- 
ditional peaks observed for the superconducting sam- 
ples BR 16 and BR 26, suggesting a second crystallo- 
graphic phase coexisting with the main phase. The co- 



existence of different crystallographic phases, was ear- 
lier reported for non-superconducting Ki_a;Fe2-ySe2i^ 
However, the orthorhombic symmetry Pmna of the ad- 
ditional phase in Ref. [SlJ is incompatible with the sym- 
metry of the second crystallographic phase observed in 
our SC samples. The phase separation in intercalated 
iron chalcogenides is now actively investigated. Recent 
transmission electron microscopy^SS x-ray diffraction,— 
ARPES,'^'* optical spectroscopy,'^ and scanning tunnel- 
ing microscopy studies^i^ provided evidence for the 
phase separation in K-intercalated iron superconduc- 
tors. Very recently, inelastic neutron scattering,'^^ optical 
spectroscopy ,-2S and Mossbauer studies^ of the supercon- 
ducting Rbo.8Fe1.6Se2 samples (batch BR 16) confirmed 
the phase separation scenario also for the Rb-based su- 
perconductors. However, the relevance to phase separa- 
tion of the second phase observed in the x-ray spectra 
of the superconducting samples is still an open question 
and additional studies are necessary to clarify its origin. 

Since the quality of the refinement of the x-ray powder- 
diffraction data within the I4/m symmetry was lower 
compared to the I4/mmm symmetry due to reduced 
statistics and a larger number of refined parameters, we 
present the data of the refinement within the I4/mmm 
symmetry. The analysis of the obtained data shows 
only slight variations of the lattice parameter a with 
the Fe concentration in the range 1.53-1.6 and a de- 
creasing tendency below 1.5 and above 1.6. An opposite 
trend is found for the variation of the c-parameter. The 
largest value of the lattice parameter a and the largest 
value of the unit cell volume are observed for the non- 
superconducting sample BR 19 with Fe concentration 
~ 1.5. With exception of this sample, the values of the 
c-parameter are lower for the superconducting samples, 
compared to the non-superconducting ones. Surprisingly, 



within the measured range of variation of the Rb concen- 
tration in the samples (0.6-0.8) no correlation between 
the Rb concentration and the value of the c-parameter 
was found, although one would intuitively expect an in- 
crease of the c-parameter with increasing Rb content. 
This probably can be attributed to a difference in the de- 
fect structure of the samples, which requires a more com- 
prehensive study. To solve remaining issues concerning 
the defect structure, vacancy ordering, and the symme- 
try of the additional crystallographic phase, single crystal 
x-ray investigations are desirable. 

Finally, we are also documenting the results of the mi- 
crostructural study. We observed that samples with the 
Fe concentration in the range of 1.53-1.6 have very flat 
surfaces without visible defects or precipitates. They can 
be easily cleaved along the a6-plane reflecting the two- 
dimensional structure of the crystals. The samples with 
the Fe content below 1.45 are also rather easily cleaved; 
their surface contains a large number of dot-like defects. 
The microstructure of the samples with the Fe concen- 
tration above 1.6 is considerably different compared to 
that of the samples with a lower Fe content. They are 
mostly three-dimensional and therefore difficult to cleave 
along the a&-plane. 



C. Magnetic susceptibility 

Figures 3 and 4, respectively, present the magnetic sus- 
ceptibility versus temperature for superconducting and 
non-superconducting samples from different growth runs 
measured in a magnetic field H=10 kOe applied paral- 
lel to the c-axis. The measurements were performed on 
cooling using a SQUID magnetometer MPMS 5 (Quan- 
tum Design). On lowering temperature from 400 K to 
33 K the susceptibility for all SC samples continuously 
decreases. The SC samples prepared by self-flux from 
the starting composition Rbo.8Fe2Se2 exhibit a higher 
susceptibility than those prepared by Bridgman method, 
which probably can be attributed to the presence of a 
small amount of non-reacted iron, e.g., of ~ 0.3% esti- 
mated from the magnetic moment for the sample F 266 
measured at 40 K that exhibits saturation at low fields. 
In contrast, the magnetization of the SC samples grown 
by Bridgman method shows non-saturated behavior (see 
below) and the amount of the non-reacted iron in these 
samples appears to be significantly reduced. In the SC 
samples a clear downturn of x at the transition temper- 
ature at around 30 K is evidenced, while x for the non- 
superconducting samples shows a small upturn towards 
lower temperatures. It is important to note the very 
close values and essentially similar temperature depen- 
dence of the susceptibility at temperatures above 33 K 
both for the SC samples from batch BR 16 and non SC 
samples from batch BR 19 despite a clear difference in 
their chemical composition. The susceptibility of the ex- 
ceptional sample BR 17 (with the strongly reduced Fe 
concentration below 1.46) shows a nonmonotonic tem- 



perature dependence with a broad maximum at around 
T*=233 K (for H\\c), which is also evidenced in the spe- 
cific heat (see below). 

In Fig. 5, the temperature dependences of the mag- 
netic susceptibility are shown for two selected SC and 
non SC samples (BR 16 and BR 19) measured in a mag- 
netic field _ff=10 kOe applied parallel and perpendicular 
to the c-axis. The susceptibility xi- for the perpendicular 
configuration (77_Lc) is by a factor of five higher than for 
H\\c, and shows a weak temperature dependence unlike 
the significant decrease of x\\ on decreasing temperature. 
Such behavior of the x± ^^'^ X\\ is characteristic for an 
anisotropic antiferromagnet, with the c-axis being the di- 
rection of alignment of the spins. A similar anisotropy of 
the magnetic properties has been recently reported for in- 
sulating TlFei.6Se2^ with a Neel temperature of 430 K. 
The anisotropic antiferromagnetism appears to be a dis- 
tinct feature of the whole family of the intercalated iron 
chalcogenides. 

In order to determine the magnetic ordering tem- 
perature we measured the susceptibility of the sample 
BR 16 at temperatures up to 650 K using a heating op- 
tion of the SQUID. The measurements (not shown) re- 
vealed a poorly pronounced maximum in the susceptibil- 
ity at a temperature of 527(±1) K. At the same time, a 
strongly irreversible behavior of x{T) on successive cool- 
ing and warming cycles was detected. Subsequent low- 
temperature measurements of the sample after the high- 
temperature cycling showed a strongly reduced amount 
of the SC phase indicating a partial degradation of the 
sample. Therefore it is unclear whether 527 K represents 
a true Neel temperature for this sample. Thus, to get an 
independent estimate of the Neel temperature we used an 
extrapolation of the x{T) data for the two measured con- 
figurations assuming that at Tj^ they should merge. The 
data on x± and xii above 100 K were fitted by polynomi- 
als of second order. From the extrapolation of x± and X|| 
a Neel temperature of about 620 (±5) K was estimated 
for SC sample BR 16. 

In Figs. 6 and 7, respectively, the magnetization curves 
M=i{H) for the superconducting sample BR 16 and non- 
superconducting sample BR 19, measured at different 
temperatures are presented. The data are shown for the 
magnetic field applied parallel and perpendicular to the 
c-axis. In the SC sample, M|| shows a small non-linearity 
at low fields and a linear increase at high fields. M|| in- 
creases significantly with temperature from 40 to 400 K. 
At the same time, M± being considerably higher than 
M|[, shows a weak temperature dependence in the tem- 
perature range 40-400 K in sharp contrast to M||. It is 
worth noting here that the magnetization of the non SC 
sample exhibits a field and temperature dependence for 
both field configurations very similar to that observed 
for the SC sample. Importantly, the magnetization of 
the non SC sample is not changing in the range from 40 
to 2 K as demonstrated by the coinciding data for 2 and 
40 K in Fig. 7. A close similarity of the magnetization 
behavior observed for both the SC and non SC samples 
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FIG. 3. (color online) Temperature dependences of the field- 
cooled susceptibility for superconducting Rbi_a;Fe2-i/Se2 
samples from different batches measured in a field of 10 kOe 
applied along the c-axis. The arrow indicates the SC transi- 
tion temperature. 
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FIG. 4. (color online) Temperature dependences of the field- 
cooled susceptibility for non-superconducting Rbi_a;Fe2-i/Se2 
samples from different batches measured in a field of 10 kOe 
applied along the c-axis. The arrow marks the temperature 
of the maximum of the susceptibility T* for sample BR 17. 



suggests that the static magnetic order coexists with the 
superconducting state. 

Figure 8 shows the temperature dependences of the 
zero-field cooled (ZFC) and field-cooled (FC) susceptibil- 
ity X for superconducting samples from different batches 
measured in a field H—IQ Oe applied parallel to the c- 
axis. The sharpest transition (ATc = 1-5 K, estimated 
from FC data) and the highest temperature of the on- 
set of superconductivity, r°"= 32.4 K was found for the 
samples from batches F 295, BR 16, and BR 26. Al- 
though the EPMA data presented above indicate some 
notable variations in the Fe and Rb concentrations for 
the samples from these batches, no significant broaden- 
ing of the SC transition or variation of T°" was detected. 
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FIG. 5. (color online) Temperature dependences of the field- 
cooled susceptibility for two SC (BR 16) and non SC (BR 19) 
Rbi-_iFe2-t/Se2 samples measured in a field of 10 kOe applied 
parallel and perpendicular to the c-axis. 
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FIG. 6. (color online) Magnetization curves for the SC sample 
BR 16 for different temperatures measured with the magnetic 
field applied along the c-axis and within the a6-plane. 



At the same time, for the SC samples from the other 
batches (for example, F 266 and BR 18), the value of 
T°" was by ~2 K lower and the variation of Xzfc indi- 
cates a broader transition compare to the samples from 
batches F 295, BR 16, and BR 26. The reduction of the 
transition temperature and broadening of the transition 
for these samples probably can be attributed to disor- 
der effects caused, for example, by larger deviations of 
the Rb and, especially, of the Fe concentrations from the 
"optimal" values of 0.8 and 1.6, respectively. 

These "optimal" values of Rb and Fe concentrations 
correspond to a formal oxidation state for the Fe ions of 
-1-2.0 as found in other superconducting Fe chalcogenides. 
Within the experimentally determined ratio of Rb to Fe 
concentration (see Table 1) we however did not found a 
change from the hole to electron doping with the varia- 
tion of the stoichiometry. Indeed, for the samples with 
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FIG. 7. (color online) Magnetization curves for the non SC 
sample BR 19 measured at different temperatures with the 
magnetic field applied along the c-axis and within the ah- 
plane. Data for 2 K (full squares) and 40 K (open squares) 
are almost coinciding. 



FIG. 8. (color online) Temperature dependences of 

ZFC and FC susceptibilities for different superconducting 
Rbi_iFe2-ySe2 samples measured in a field of 10 Oe applied 
along the c-axis. The arrow indicates the temperature of the 
onset of the superconducting transition r°" = 32.4 K for 
sample BR 16. 



the lowest Fe concentration (below 1.5, e.g., BR17) the 
calculated valence of Fe is -(-2.27. For the samples with 
high Fe content (above 1.61, e.g., BR22) the valence of Fe 
is +2.07. Measurements of the thermoelectric power are 
necessary to clarify the role of doping on the properties 
of the Rb-based iron chalcogenides, as it was shown for 
the related K-intercalated compounds 42ii^ 

It is important to note that the values of xfc (Meiss- 
ner effect) at 2 K for the SC samples are at least three 
order of value lower than the values of xzfc (shielding 
effect), probably due to strong pinning. At the same 
time, the values of 47rx at 2 K calculated from the xzfc 
data are far above unity pointing to the dominance of 
demagnetizing effects. The measurements of xzfc for 
the in-plane configuration (with R perpendicular to the 
c-axis) with the negligible demagnetizing factor for the 
samples from batches F 295, BR 16, and BR 26 gave 
a value of 47rx close to unity that may indicate a large 
volume fraction of the SC phase which yields a perfect 
100 % shielding. 

Figure 9 presents the magnetization hysteresis loops 
for SC samples measured at 2 K with the magnetic field 
H applied along the c-axis. The hysteresis loop for this 
configuration is symmetric with respect to the origin 
(0.0). The diamagnetic response for samples from differ- 
ent batches differs by a factor of 20. The reason of such 
a strong difference is unclear. It can be caused by varia- 
tions of the stoichiometry detected by EPMA (see Table 
m or by changes in the preparation conditions, such as 
oxygen contamination. The phase separation can also be 
critical for obtaining high SC hysteresis width. Indeed, 
the highest diamagnetic response was found in the sam- 
ples from batches F 295, BR 16 and BR 26 in which x-ray 
diffraction revealed strong reflexes of the additional crys- 
tallographic phase. Additional pinning centers created 
by this additional crystallographic phase can enhance the 



diamagnetic response. 

The critical current density j^, at 2 K estimated 
from the hysteresis data for the samples with the high- 
est SC parameters using the Bean model for hard 
superconductors^^i^ reaches a value of 1.6 x 10''' Kjcvc? . 
It is by a factor of five lower than the intrinsic value of 
jc reported for FeSe substituted with Te4^i^ 

In Fig. 10 the hysteresis loop for the SC sample BR 16 
measured at 2 K with the magnetic field applied perpen- 
dicular to the c-axis is plotted. The maximum values of 
magnetization for _ff _L c and -ff || c correlate with the 
difference in the demagnetization factors for these two 
configurations. The hysteresis loop for the configuration 
if II c reveals two contributions, one superconducting, 
and the other linear in magnetic field. For comparison, in 
Fig. 10 the magnetization curve M||=f(_ff) for the same 
sample measured at 40 K, just above the SC transition, 
is shown. It agrees well with the linear contribution that 
can be derived from the data at 2 K. This again points 
out that the magnetic correlations associated with the 
static antiferromagnetic order coexist with the SC state. 



D. Resistivity 

Figures 11(a) and (b) demonstrate the temperature 
dependences of the resistivity for the SC and non SC 
samples, respectively. The resistivity measurements were 
done on rectangular samples by a four-point method in 
the temperature range 2-300 K using the resistivity mea- 
surement option of the Quantum Design physical proper- 
ties measurement system (PPMS) in magnetic fields up 
to 90 kOe. Both, the in-plane resistivity, pah-, with cur- 
rent i in the a&-plane, and the inter-plane resistivity, pc, 
with current j parallel to the c-axis were measured. The 
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FIG. 9. (color online) Hysteresis loops measured at 2 K with 
magnetic field applied along the c-axis for different SC sam- 
ples Rbi_iFe2-ySe2. For clarity, the data for samples F 266 
and BR 18 are magnified by factors of 5 and 10, respectively. 




FIG. 10. (color online) Hysteresis loop for SC sample BR 16 
at 2 K (solid circles) measured with the magnetic field applied 
perpendicular to the c-axis, together with the magnetization 
curve at 40 K (open squares). 



magnetic field was always perpendicular to the current 
direction, so only the transverse configuration was stud- 
ied. The contacts were made with a conductive silver 
paint. The error in the absolute value of the measured 
resistivity was about 20 % due to finite dimensions of 
the potential contacts. The in-plane resistivity pat of the 
SC sample prepared by self-flux method (F 266) is by 
approximately one order of magnitude larger than that 
for samples grown by Bridgman method which can be at- 
tributed to higher impurity content of the self-flux sam- 
ples. The inter-plane resistivity pc of a typical sample 
grown by Bridgman method with high SC parameters 
(BR 16) is by a factor of two higher than the in-plane 
resistivity. The resistivity curves of all superconducting 
samples show a broad hump at a temperature T^ax that 
for different samples varies in the range of 190-215 K 
with a semiconductor-like temperature dependence above 



and a metal-like dependence below Tmax- We note that 
the superconducting samples from batch BR 16 have 
the smallest value of the resistivity compared to other 
batches. It is also the smallest for the superconduct- 
ing Rbi_j:Fe2-ySe2 crystals reported so far.^^'^^ Together 
with the larger resistivity ratio p{Tmax) I p{Tc) = 37 this 
may indicate a higher purity of these superconducting 
samples. 
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FIG. 11. (color online) (a) Temperature dependences of the 
in-plane resistivity pab and inter-plane resistivity pa for the 
SC sample BR 16 (left scale) and for the in-plane resistivity 
for the SC sample F 266 (right scale) measured on cooling in 
zero external magnetic field, (b) Temperature dependences 
of the in-plane resistivity for the non SC samples. Dashed 
curves mark the Arrhenius fits at high temperatures, solid 
lines - the MVRH fits al lower temperatures as described in 
the text. 

The resistivity of the non-superconducting samples is 
by several orders higher than that of the SC samples. 
It shows a steep increase on decreasing temperature and 
a small but clearly discernible anomaly at T*=233 K. 
In the range 300-240 K the temperature dependence of 
the resistivity for the non SC samples can be described 
by Arrhenius law p = pocxp{AEa/kT) with an activa- 
tion energy AE^ of 0.16, 0.079 and 0.073 eV for the 
samples BR 19, BR 22, and BR 17, respectively. At 
temperatures below 240 K down to 80 K the resistiv- 
ity of the non SC samples is not thermally activated 



but, instead, can be reasonably described within a Mott 
variable-range hopping (MVRH) model by an expres- 
sion p = poexp(To/T)^/'*. The resistivity of the sample 
BR 22 with the Fe concentration above 1.62 is also show- 
ing semiconductor-like behavior with decreasing temper- 
ature, but its value is by several orders of magnitude 
lower than that of the samples with the Fe concentration 
below 1.5. Summarizing the observed resistivity behavior 
of the non-superconducting samples we notice rather low 
values of the activation energies which are more typical 
for heavily doped semiconductors or Mott insulators than 
for intrinsic band insulators. Additional studies including 
Hall effect are necessary to clarify the actual transport 
mechanisms in the non-superconducting Rbi_a;Fe2_j,Se2 
crystals. 

Figures 12(a) and (b), respectively, illustrate the ef- 
fect of a transverse magnetic field on the in-plane pab 
and inter-plane pc resistivities in the transition region 
for the SC sample BR 16. The measurements were done 
on warming after cooling in zero field. In zero field the 
transition temperature, determined at the level of the 
90 % drop of the normal state resistivity /5„, coincides 
with the r°" (within 0.1 K) determined from the sus- 
ceptibility measurements. This again indicates the high 
quality of the superconducting samples. When increas- 
ing the magnetic field, the resistivity curves are displaced 
to lower temperatures. The inset in Fig. 12(b) shows 
the temperature dependence of the upper critical field 
Hc2{T) estimated using the criterion of the 90 % drop 
of the normal state resistivity. On approaching Tc, the 
slope of the if c2 (T) curve becomes smaller compared to 
that at lower temperatures. The estimation of the up- 
per critical field i?c2(0) for T = K made within the 
Werthamer-Helfand-Hohenberg model^ using the rela- 
tion Hc2{0) = -0.69TeidHc2{T)/dT) |t, gave a value of 
250 kOe for the in-plane and 630 kOe for the inter-plane 
configurations. They are by a factor of two lower than 
those for FeSe substituted with Te^i^ and are the lowest 
reported so far for SC Rbi_2;Fe2-ySe2 samplesJ^ii^ 



E. Specific heat 

Figure 13 presents the temperature dependences of the 
specific heat C for selected superconducting and non- 
superconducting samples. The heat capacity was mea- 
sured by a relaxation method using the PPMS in the 
temperature range 1.8-300 K and in magnetic fields up 
to 90 kOe. The magnetic field was applied parallel to 
the c-axis of the samples. The specific heat for different 
samples exhibits very similar behavior and values in the 
measured temperature range. The specific heat of the 
non SC sample BR 17 shows a sharp peak at around 
220 K. Its position is lower than that of the peak in 
the susceptibility at 233 K for the same field configu- 
ration. Interestingly, the susceptibility of this sample 
for the field configuration perpendicular to the c-axis ex- 
hibits an anomaly exactly at the same temperature of 
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FIG. 12. (color online) (a) Temperature dependences of the 
in-plane resistivity pab for the sample BR 16 measured in var- 
ious magnetic fields applied parallel to the c-axis. (b) Tem- 
perature dependences of the inter-plane resistivity pc for the 
sample BR 16 measured in various magnetic fields applied 
perpendicular to the c-axis. Inset: temperature dependence 
of the upper critical field for inter-plane (closed symbols) and 
in-plane (open symbols) configurations as described in the 
text. 



220 K as observed in the specific heat. The application 
of the magnetic field of 50 kOe does not have any in- 
fluence on the temperature position of the peak in the 
specific heat at 220 K indicating that it is not related to 
conventional AFM or FM transitions. Temperature de- 
pendent x-ray measurements in the range 300-100 K did 
not find any evidence of a structural transformation ei- 
ther. Preliminary Mossbauer studies of the samples from 
this batch revealed new lines below T* typical for triva- 
lent Fe ions4S The observed anisotropy in the value of 
T* for i7|| c-axis and _ff_Lc-axis suggests a very peculiar 
type of this transition which needs further investigation. 

In Fig. 14 the temperature dependences of the spe- 
cific heat is shown in the representation C/T vs. T^ for 
different samples at temperatures below 12 K. A linear 
dependence with two different slopes below and above 
5 K is clearly discernible. The values of the residual 
Sommerfeld coefficient jr were determined from the fit 
to the experimental data at temperatures below 3 K us- 
ing the expression C/T = jj. + PT'^ , where the prefactor 
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FIG. 13. (color online) Temperature dependences ol the spe- 
cific heat for different Rbi-a;Fe2-j/Se2 samples. The arrow 
marks the anomaly in the specific heat at T*=220 K for the 
sample BR 17. 



/3 characterizes both the lattice and magnon contribu- 
tions to the specific heat, which are proportional to T^, 
and therefore cannot be separated. The calculated val- 
ues of these parameters are given in Table |TT1 We notice 
a rather low value of jr obtained for the superconduct- 
ing samples, which indicates their high purity and mini- 
mal structural disorder or impurity content. The values 
of 7r vary from 0.18 to 0.85 mJ/molK^ for the samples 
with high SC (BR 16) and poor SC parameters (F 266), 
respectively, reflecting higher impurity content and dis- 
order of the poorly SC samples in agreement with the 
susceptibility results. Similar values of 7,- were reported 
for the related SC Ki_2;Fe2-j/Se2 compound. ^° At the 
same time, the non-superconducting samples BR 17 and 
BR 22 exhibit rather high values of 7^, which probably 
can be attributed to structural disorder due to Fe and Rb 
vacancies. Importantly, the non-superconducting sam- 
ple BR 19 exhibits an extremely low value of 7^ indicat- 
ing minimal disorder, being correlated with its best crys- 
tallinity revealed by x-ray studies (see Fig. 2). It should 
be noted that the value of f3 for the sample BR 19 is 
the closest to those of the superconducting samples. To- 
gether with very similar values of the susceptibility and 
its temperature dependence (see Fig. 5), this justifies the 
use of the specific heat data of the sample BR 19 as a 
reference for the lattice and magnon contributions. The 
respective contributions for the superconducting samples 
were calculated using the specific heat data of the sam- 
ple BR 19 scaled by the square root of the ratio of their 
molar mass. 

Fig. 15 shows the electronic specific heat in the repre- 
sentation Ce/T VS. T for several superconducting sam- 
ples from different batches obtained by subtraction of 
the lattice and magnon contributions from the total mea- 
sured specific heat. It has been found that the supercon- 
ducting contribution to the specific heat is rather small 
compared to that of the lattice and magnon contribu- 
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14. (color online) Temperature dependences of the spe- 

in the representation C jT vs. T^ for different 

-j,Se2 samples at temperatures below 12 K. Upper 

T in different applied magnetic fields for the 

Lower inset: magnetic field dependence of the 
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residual Sommerfeld coefficient 7r for different superconduct- 
ing samples with high SC parameters. 



tions. However, a pronounced lambda-like anomaly of Ce 
at the superconducting transition is clearly evidenced for 
all SC samples. The sample F 266 with poor SC param- 
eters exhibits a broad anomaly of Ce at T^-, which is dis- 
placed to lower temperatures compared to the anomaly 
in samples from the batch BR 16 with high SC parame- 
ters. The samples from batch BR 26 show the sharpest 
anomaly of Ce at T^.. 
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FIG. 15. (color online) Temperature dependences of the elec- 
tronic specific heat Ce/T vs. T for different Rbi_a;Fe2-i/Se2 
samples calculated by subtraction of the magnon and lattice 
contributions. The arrows mark Tc for the samples F 266 and 
BR 16. 

The value of the Sommerfeld coefficient in the normal 
state 7„ for the SC samples was estimated from the con- 
straint of entropy conservation at the onset of Tc, i.e.. 



pTe 



/p " Ce/TdT = Jq" JndT. The values of 7„ are given 
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TABLE II. Parameters determined from the susceptibility {T°", ATc) and heat capacity measurements. 



Sample 


rpon 

(K) 


ATc 

(K) 


7r 

(mJ/mol K^) 


P 
(mJ/mol K-*) 


7n 

(mJ/mol K^) 


AC/7„r, 


F266 


30.9 


2.8 


0.85 


0.96 


10 


1.0 


BR 16 


32.4 


1.5 


0.18 


1.09 


12 


0.95 


BR 16-8 


32.4 


1.6 


0.46 


1.04 


12 


0.96 


BR 17 


- 


- 


0.33 


0.90 


- 


- 


BR 19 


- 


- 


- 


0.95 


- 


- 


BR 22 


- 


- 


3.0 


1.36 


- 


- 


BR 26 


32.4 


1.5 


0.34 


0.94 


7.4 


2.2 



in Table II. They vary between 7.4 and 12 mJ/molK^ 
for SC samples from different batches and are compara- 
ble to those reported for SC Ki_a;Fe2-7/Se2i^ Assuming 
that the residual Sommerfeld coefficient 7^ corresponds 
to the fraction of the normal conducting state, the ob- 
tained ratio of ^rlln implies that the volume fraction of 
the superconducting phase in our samples varies within 
92-98 %. The values of the jump of the specific heat 
AC/Tc at Tc were calculated by subtracting the value of 
7„ from the maximum value of Ce/T at Tc- They vary 
from 10 to 16 mJ/molK^ for different SC samples and are 
by a factor of '-^ 5 smaller than the respective values of 
AC/Tc observed for single crystaUine FeSeo.sTeo.s^l'^'^ 
The reduced specific heat jump at Tc, IS.C /^nTc varies 
from 0.95 (for BR 16) to 2.25 (for BR 26). The values 
below 1.0 are probably related to errors in the estimation 
of 7„ for these samples. The obtained value of AC/7„Tc 
= 2.25 for the sample BR 26 with the sharpest anomaly in 
the specific heat at Tc is higher than the BCS estimate of 
1.43 for the weak-coupling limit and may indicate strong- 
coupling superconductivity in Rbi_a;Fe2-j,Se2 similar to 
that observed in Te substituted FeSe^ and in optimally 
doped Bai_:rK^Fe2As2.^"-^ 

Since specific heat in the SC regime varies substantially 
among samples with rather close stoichiometry, we did 
not attempt to estimate the value of the superconduct- 
ing gap. However, to get an insight into the symmetry 
of the order parameter we measured the low-temperature 
specific heat in different magnetic fields. The upper in- 
set in Fig. 14, where the specific heat in the repre- 
sentation C /T vs. T^ is plotted for one of the sam- 
ples with high SC parameters (BR 16), demonstrates 
the effect of the magnetic field on the specific heat at 
lowest temperatures. The application of the magnetic 
field results in an increase of the residual Sommerfeld 
coefficient 7^ probably related to the pair-breaking ef- 
fect. The field dependences of 7^ for three samples 
with high SC parameters are shown in the lower inset 
in Fig. 14. The values of 7^ for these samples are very 
close and exhibit a similar field dependence that signif- 
icantly differs from the linear dependence expected for 
fully gapped superconductors. It must be noted that the 
dependence of ^r{H) is close to \/H characteristic for 
(i-wave symmetry.— Recently, a Volovik-like ^/H term 



in the specific heat was reported for an optimally doped 
BaFe2(Aso.7Po.3)2 pnictide which indicate the presence 
of nodes in the superconducting gapiS However, there 
is no experimental evidence of nodes in Ai_2;Fe2-2/Se2 
yet. On the contrary, very recent ARPES studies of re- 
lated Ai_2;Fe2-j/Se2 materials22i2^ seem to fully rule out 
d-wave pairing. At present, the sub-linear dependence 
of the jr{H) observed in Rbi_a;Fe2-2/Se2 is lacking any 
reasonable theoretical explanation but multiband effect 
with two different gaps discussed in Ref. [53 can be also 
relevant in this case. 



III. SUMMARY AND CONCLUSIONS 

Our detailed preparative, structural, magnetic, con- 
ductivity, and specific heat studies of Rbi_a:Fe2-i,Se2 sin- 
gle crystals testify several important peculiarities of these 
materials; 

1. Three different regions with distinct structural, 
magnetic, and conducting behavior of the samples are 
documented in the Rbi_a;Fe2-ySe2 system as presented 
in the phase diagram in Fig. 16. Superconductivity is ob- 
served in the range of the Fe concentration above 1.53 and 
below 1.60 (region II of the diagram). Superconducting 
samples with the sharpest transitions, the highest crit- 
ical temperature of 32.4 K and the largest diamagnetic 
response have the "optimal" composition Rbo.8Fe1.6Se2. 
The reduction of the superconducting parameters of the 
samples can be attributed to deviations of the Fe con- 
centration from the "optimal" value of 1.6. Within the 
experimentally determined range of the Rb concentration 
(0.6-0.8) no correlation between the variation of the Rb 
content and the lattice parameters of the samples was 
found. The superconducting samples show the smallest 
value of the lattice parameter c compared to the non- 
superconducting samples. Samples with the Fe content 
below 1.5 (region I) and above 1.6 (region HI) show an 
insulating, or semiconducting behavior at low tempera- 
tures, respectively. 

2. In the whole range of the studied Fe concentra- 
tions the susceptibility measurements evidenced a be- 
havior typical for an anisotropic antiferromagnet with a 
transition temperature T/v above 500 K. 
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FIG. 16. (color online) Phase diagram of Rbi_2;Fe2_ySe2 sys- 
tem presenting dependences of the lattice constant c, transi- 
tion temperature Tc and the Neel temperature Tn for three 
different regions with distinct structural, magnetic and con- 
ducting behavior. 



3. Magnetization measurements provide an evidence 
for the presence of antiferromagnetic correlations in the 
ordered superconducting state in samples with the Fe 
concentration in the range of 1.53-1.6. This suggests the 
coexistence of superconductivity and static antiferromag- 
netic order. In a view of recent experiments^"— this 
coexistence should be considered as being on the macro- 
scopic level within a scenario of the phase separation. 

4. Magnetic hysteresis measurements of the supercon- 
ducting samples revealed high values of the critical cur- 
rent density jc of 1.6 x 10'* A/cm^ for superconducting 
samples with the composition close to Rb2Fe4Se5. The 
obtained value of jc is the largest reported thus far for 
Ai_j;Fe2_2/Se2 materials. 



5. The upper critical fields Hc2 of ~ 250 kOe for the 
in-plane, and 633 kOe for the inter-plane configurations 
are estimated from the resistivity study of the supercon- 
ducting samples in magnetic fields. They are the lowest 
reported to date for the Rbi_a;Fe2_j,Se2 superconductors. 

6. Specific-heat measurements of the superconduct- 
ing samples evidenced very low values of the residual 
Sommerfeld coefficient corresponding to 92-98% volume 
fraction of the superconducting phase, which confirm the 
bulk nature of the superconducting state. 

7. Resistivity studies revealed an insulating behavior 
at low temperatures for samples with Fe concentrations 
below 1.5. The conductivity of these samples at low tem- 
peratures can be reasonably described by the variable- 
range hopping Mott mechanism. The conductivity of the 
samples with Fe concentrations above 1.6 is typical for 
heavily doped semiconductors. 

8. In the non-superconducting samples with the Fe 
concentration below 1.45 an unusual phase transition be- 
low 233 K is revealed which is not related to antiferro- 
magnetic or structural transformations. 

9. The discrepancies in the range of Fe concentrations, 
where superconductivity shows up in Rbi_a;Fe2-ySe2 as 
compared to earlier reports, are likely to be related to 
the accuracy of the methods used for the compositional 
analysis. 



ACKNOWLEDGMENTS 

The authors thank Dana Vieweg for experimental sup- 
port. This research has been supported by the DFG 
via SPP 1458 under Grant DE1762/1-1 and the Transre- 
gional Collaborative Research Center TRR 80 (Augsburg 
- Munich). 



* Y. Kamihara, T. Watanabe, M. Hirano, H. Hosono, J. Am. 

Chem. Soc. 130, 3296 (2008). 
^ M. Rotter, M. Tegel, D. Johrendt. Phys. Rev. Lett. 101, 

107006 (2008). 
^ J.H. Tapp, Z. Tang, B. Lv, K. Sasmal, B. Lorenz, C.W. 

Chu, A.M. Guloy, Phys. Rev. B 78, 060505 (2008). 
^ F.-C. Hsu, J.Y. Luo, K.W. Yeh, T.K. Chen, T.W. Huang, 

P.M. Wu, Y.-C. Lee, Y.L. Huang, Y.-Y. Chu, D.C. Yan, 

M.-K. Wu, Proc. Natl. Acad. Sci. U.S.A. 105, 14262 

(2008). 
^ S. Medvedev, T.M. McQueen, LA. Troyan, T. Palasyuk, 

M.I. Eremets, R.J. Cava, S. Naghavi, F. Casper, V. Kseno- 

fontov, G. Wortmann, and C. Felser, Nature Materials 8, 

630 (2009). 
® S. Margadonna, Y. Takabayashi, Y. Ohishi, Y. Mizuguchi, 

Y. Takano, T. Kagayama, T. Nakagawa, M. Takata, and 

K. Prassides, Phys. Rev. B 80, 064506 (2009). 
^ M.H. Fang, H.M. Pham, B. Qian, T.J. Liu, E.K. Vehstedt, 

Y. Liu, L. Spinu, and Z.Q. Mao, Phys. Rev. B 78, 224503 

(2008). 



* K.W. Yeh, T.-W. Huang, Y. Huang, T.-K. Chen, F.-C. 

Hsu, P.M. Wu, Y.-C. Lee, Y.-Y. Chu, C.-L. Chen, J.-Y. 

Luo, D.-C. Yan, and M.-K. Wu, Europhys. Lett. 84, 37002 

(2008). 
^ J. Guo, S. Jin, G. Wang, S. Wang, K. Zhu, T. Zhou, M. 

He, and X. Chen, Phys. Rev. B 82, 180520 (2010). 
*° A. Krzton-Maziopa, Z. Shermadini, E. Pomjakushina, V. 

Pomjakushin, M. Bendele, A. Amato, R. Khasanov, H. 

Luetkens, and K. Conder, J. Phys. Condensed Matter 23, 

052203 (2011). 
" J. J. Ying, X. F. Wang, X. G. Luo, A. F. Wang, M. Zhang, 

Y. J. Yan, Z. J. Xiang, R. H. Liu, P. Cheng, G. J. Ye, and 

X. H. Chen, Phys. Rev. B 83, 212502 (2011). 
*2 C.H. Li, B. Shen, F. Han, X. Zhu, and H.H. Wen, Phys. 

Rev. B 83, 184521 (2011). 
*^ X. G. Luo, X. F. Wang, J. J. Ying, Y. J. Yan, Z. Y. Li, 

M. Zhang, A. F. Wang, P. Cheng, Z. J. Xiang, G. J. Ye, 

R. H. Liu, X. H. Chen, New J. Phys. 13, 053011 (2011). 
" Y. Mizuguchi, H. Takeya, Y. Kawasaki, T. Ozaki, S. 

Tsuda, T. Yamaguchi, Y. Takano Appl. Phys. Lett. 98, 



12 



042511 (2011). 
^^ A. F. Wang, J. J. Ying, Y. J. Yan, R. H. Liu, X. G. Luo, 

Z. Y. Li, X. F. Wang, M. Zhang, G. J. Ye, P. Cheng, Z. J. 

Xiang, and X. H. Chen, Phys. Rev. B 83, 060512 (2011). 
^^ H.D. Wang, C.H. Dong, Z.J. Li, Q.H. Mao, S.S. Zhu, CM. 

Feng, H. Q. Yuan, and M.H. Fang, Europhys. Lett. 93, 

47004 (2011). 
^"^ Z. Shermadini, A. Krzton-Maziopa, M. Bendele, R. 

Khasanov, H. Luetkens, K. Gender, E. Pomjakushina, S. 

Weyeneth, V. Pomjakushin, O. Bossen, and A. Amato, 

Phys. Rev. Lett. 106, 117602 (2011). 
^^ L. Li, Z. R. Yang, Z. T. Zhang, W. Tong, C. J. Zhang, S. 

Tan, Y. H. Zhang, arXiv: 1101.5327 (unpublished). 
1^ R. H. Liu, X. G. Luo, M. Zhang, A. F. Wang, J. J. Ying, 

X. F. Wang, Y. J. Yan, Z. J. Xiang, P. Cheng, G. J. Ye, Z. 

Y. Li and X. H. Chen, Europhys. Lett. 94, 27008 (2011). 
2° M.H. Fang, H.D. Wang, C.H. Dong, Z.J. Li, CM. Feng, J. 

Chen, and H. Q. Yuan, Europhys. Lett. 94, 27009 (2011). 

21 W. Bao, Q. Huang, G. F. Chen, M. A. Green, D. M. Wang, 
J. B. He, X. Q. Wang, Y. Qiu, arXiv; 1102.0830 (unpub- 
lished) . 

22 D. Mou, S. Liu, X. Jia, J. He, Y. Peng, L. Zhao, L. Yu, 
G. Liu, S. He, X. Dong, J. Zhang, H. Wang, C Dong, M. 
Fang, X. Wang, Q. Peng, Z. Wang, S. Zhang, F. Yang, 
Z. Xu, C Chen, and X. J. Zhou, Phys. Rev. Lett. 106, 
107001 (2011). 

2=* L. Zhao, D. Mou, S. Liu, X. Jia, J. He, Y. Peng, L. Yu, 
X. Liu, G. Liu, S. He, X. Dong, J. Zhang, J. B. He, D. M. 
Wang, G. F. Chen, J. G. Guo, X. L. Chen, X. Wang, Q. 
Peng, Z. Wang, S. Zhang, F. Yang, Z. Xu, C Chen, and 
X. J. Zhou, Phys. Rev. B 83, 140508 (2011). 

24 T. Qian, X.-P. Wang, W.-C Jin, P. Zhang, P. Richard, G. 
Xu, X. Dai, Z. Fang, J.-G. Guo, X.-L. Chen, and H. Ding, 
Phys. Rev. Lett. 106, 187001 (2011). 

25 I. Mazin, Physics 4, 26 (2011). 

2"^ J. Rodriguez-Carvajal, Physica B 192, 55 (1993). 

2'' Z. Gao, Y. Qi, L. Wang, C Yao, D. Wang, X. Zhang, Y. 

Ma, arXiv: 1103.2904 (unpublished). 
2« P. Zavalij, Wei Bao, X. F. Wang, J. J. Ying, X. H. Chen, D. 

M. Wang, J. B. He, X. Q. Wang, G. F. Chen, P.-Y. Hsieh, 

Q. Huang, and M. A. Green, Phys. Rev. B 83, 132509 

(2011). 
2^ V. Yu. Pomjakushin, D. V. Sheptyakov, E. V. 

Pomjakushina, A. Krzton-Maziopa, K. Conder , D. 

Chernyshov, V. Svitlyk, and Z. Shermadini, Phys. Rev. 

B 83, 144410 (2011). 
^° F. Ye, S. Chi, Wei Bao, X. F. Wang, J. J. Ying, X. H. Chen, 

H. D. Wang, C H. Dong, Minghu Fang, arXiv: 1102.2882 

(unpublished) . 
31 W. Bao, G. N. Li, Q. Huang, G. F. Chen, J. B. He, M. A. 

Green, Y. Qiu, D. M. Wang, J. L. Luo, Phys. Rev. Lett. 

107, 137003 (2011). 
=^2 Z. Wang, Y. J. Song, H. L. Shi, Z. W. Wang, Z. Chen, H. 

F. Tian, G. F. Chen, J. G. Guo, H. X. Yang, and J. Q. Li, 

Phys. Rev. B 83, 140505 (2011). 
33 A. Ricci, N. Poccia, G. Campi, B. Joseph, G. Arrighetti, 

L. Barba, M. Reynolds, M. Burghammer, H. Takeya, Y. 

Mizuguchi, Y. Takano, M. Colapietro, N. L. Saini, and A. 

Bianconi, Phys. Rev. B 84, 060511 (2011). 



3* F. Chen,l M. Xu, Q.Q. Ge, Y. Zhang, Z.R. Ye, L.X. Yang, 

J. Jiang, B.P. Xie, R.C Che, M. Zhang, A.F. Wang, X.H. 

Chen, D.W. Shen, X.M. Xie, M.H. Jiang, J.P. Hu, and 

D.L. Feng, [ajXiv:1106.3026 (unpublished). 
35 R.H. Yuan, T. Don g, Y.J. S ong, G.F. Chen, J.P. Hu, J.Q. 

Li, and N.L. Wang, arXiv:1102.1381 (unpublished). 
3« W. Li, H. Ding, P. Deng, K. Chang, C Song, K. 

He, L. Wang, X. Ma, J.-P. Hu, X. Chen, Q.-K. Xue, 

arXiv: 1108.0069 (unpublished). 
3^ P. Cai, C Ye, W. Ruan, X. Zhou, A. Wang, M. Zhang, X. 

Chen, Y. Wang, arXiv:1108.2798 (unpublished). 
38 J. T. Park, G. Friemel, Yuan Li, J.-H. Kim, V. Tsurkan, J. 

Deisenhofer, H.-A. Krug von Nidda, A. Loidl, A. Ivanov, 

B. Keimer, D. S. Inosov, arXiv:1107.1703 (unpublished). 
3^ A. Charnukha, J. Deisenhofer, D. Prpper, M. Schmidt, 

Z. Wang, Y. Goncharov, A. N. Yaresko, V. Tsurkan, B. 

Keimer, A. Loidl, A. V. Boris, larXiv: 1108.56 98 (unpub- 
lished) . 
'"' V. Ksenofontov, G. Wortmann, S. Medvedev, V. Tsurkan, 

J. Deisenhofer, A. Loidl, C Felser, larXiv:1108.3006] (un- 
published) . 
■'i B. C Sales, M. A. McGuire, A. F. May, H. Cao, B. C 

Chakoumakos, and A. S. Sefat, Phys. Rev. B 83, 224510 

(2011). 
'*2 K. Wang, H. Lei, and C Petrovic, Phys. Rev. B 83, 174503 

(2011). 
■'3 Y. J. Yan, M. Zhang, A. F. Wang, J. J. Ying, Z. Y. 

Li, W. Qiu, X. G. Luo, J. Q. Li, J. Hu, X. H. Chen, 

arXiv:1104.4941 (unpublished). 
^^ CP. Bean, Phys. Rev. Lett. 8, 250 (1962). 
^^ CP. Bean, Rev. Mod. Phys. 36, 90 (1964). 
*® A. Giinther, J. Deisenhofer, Ch. Kant, H.-A. Krug von 

Nidda, V. Tsurkan, and A. Loidl, Supercond. Sci. Technol. 

24, 045009 (2011). 
*'' V. Tsurkan, J. Deisenhofer, A. Giinther, Ch. Kant, M. 

Klemm, H.-A. Krug von Nidda, F. Schrettle, and A. Loidl, 

Eur. Phys. J. B 79, 289 (2011). 
■'** N.R. Werthamer, E. Helfand, and P.C Hohenberg, Phys. 

Rev. 147, 295 (1966). 
"^^ V. Ksenofontov et al., to be published. 
5"^ B. Zeng, B. Shen, G. F. Chen, J. B. He, D. M. Wang, C 

H. Li, and H. H. Wen, Phys. Rev. B 83, 144511 (2011). 

51 A. Serafin, A. I. Coldea, A. Y. Ganin, M. J. Rosseinsky, 
K. Prassides, D. VignoUes, and A. Carrington, Phys. Rev. 
B 82, 104514 (2010). 

52 J. Hu, T. J. Liu, B. Qian, A. Rotaru, L. Spinu, and Z. Q. 
Mao, Phys. Rev. B 83, 134521 (2011). 

53 G. Mu, H. Luo, Z. Wang, L. Shan, C Ren, and H.-H. Wen, 
Phys. Rev. B 79, 174501 (2009). 

5'' P. Popovich, A. V. Boris, O. V. Dolgov, A. A. Golubov, 

D. L. Sun, C T. Lin, R. K. Kremer, and B. Keimer, Phys. 

Rev. Lett. 105, 027003 (2010). 
55 C. Kant, J. Deisenhofer, A. Giinther, F. Schrettle, A. Loidl, 

M. Rotter, and D. Johrendt, Phys. Rev. B 81, 014529 

(2010). 
5S G. E Volovik, JETP Lett. 58, 469 (1993). 
5^ Y. Wang, J.S. Kim, G. R. Stewart, P.J. Hirschfeld, S. 

Graser, S. Kasahara, T. Terashima, Y. Matsuda, T. 

Shibauchi, I. Vekhter.' arXKr:1109.0554, funpublished'). 



